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The effective action formahsm of quantum field theory is used to study the prop- 
erties of the non-relativistic interacting Bose gas. The Gaussian approximation is 
formulated by calculating the effective action to the first order of the optimized 
expansion. In the homogeneous limit the method respects the Hughenholz-Pines 
theorem, leading to the gapless spectrum both for excitations and for density fluc- 
tuations. Renormalization is carried out by adopting dimensional regularization. 
The results for critical temperature are compared with that obtained in the loop 
expansion and lattice calculations. 

I. INTRODUCTION 

A weakly interacting Bose gas was intensively studied in the period 1947-1965 as the 
theoretical example of a Bose-condensed system, described by non-relativistic QFT. Various 
approximation methods have been considered Q-^l, achieving a great success in describing 
the essential features of '^He superfluidity. However the interactions in quantum liquids are 
too strong and the weakly interacting system was in quest for a quantitative test of many- 
body methods. Such a system become accessible in 1995 after the experimental verification of 
the Bose- Einstein condensation (BEC) in dilute gases in magnetic traps 0] . When adapting 
many-body approximations to the trapped gases, the problem of approximation consistency 
0] attracted a renewed attention The approximation scheme should preserve relations 
between physical observables, arising from the symmetry of the theory (Ward identities). 
For spatially uniform systems, such a relation is implied by the Hughenholz-Pines (H-P) 
theorem (the non-relativistic analogue of the Goldstone theorem). The H-P theorem 
shows that the single particle spectrum of a many-body system is gapless, if a global sym- 
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metry is spontaneously broken. This imposes strong constraints on approximations. The 
simplest Bogoliubov approximation Q| fulfils the constraints, but including the first order 
corrections to propagators violates the H-P theorem. Only the second order approximation, 
developed by Beliaev 0] and extended to finite temperature by Popov, is fully consistent. 
In the self-consistent Hartree-Fock-Bogoliubov approximation the H-P is violated: there is 
a gap in the single particle spectrum although the density fluctuations are gapless 

The Bose condensed systems are usually studied within field theoretical techniques in- 
volving the particle Hamiltonian. The Lagrangian formulation is more convenient for dis- 
cussing the general issues, such as symmetries, conservation laws, and Ward identities. Upon 
straightforward extension to finite temperature the grand canonical thermodynamic poten- 
tial is obtained as a function of the superfiuid order parameter. However, the results obtained 
in the Lagrangian approach are also controversial: Toyoda P] claimed the one-loop effective 
potential at finite temperature to be consistent with the H-P theorem, but later an oppo- 
site statement has been made Q]. Methods of self-consistent re-summation have been also 
discussed but the consistency with the H-P theorem is unclear. We show that the 

effective action provides a useful tool to construct consistent approximations. 

The outline of our work is as follows. In section |H] the effective action formalism is 
introduced using path-integrals for generating functionals of non-relativistic theory. We 
show that the H-P theorem can be formulated in a form of a simple criterium, which is 
useful for checking approximations consistency. In section IIIII we study the loop expansion 
for the effective action. The one-loop result gives the Beliaev-Popov approximation at 
finite temperature. In section HVl we extend the method of the optimized expansion jl^ to 
the effective action of non-relativistic theory. The self-consistent Gaussian approximation, 
obtained in the first order, fulfils the H-P theorem. Results for the critical temperature and 
conclusions are summarized in Section IVl 



II. EFFECTIVE ACTION 

The system of many spinless atoms interacting by two-body forces is described by a 
complex scalar field $(r,t) with the Lagrangian density 
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-K.t(r)$*(r,t)$(r,t)-ly'ciV($*(r,t)$(r,t)t/(r-r')$*(r',t)$(r',t)), (1) 

where U{r — r') represents the two-body interatomic potential, and Vexti'"^) the external 
potential of the trap. The Lagrangian is invariant under the global U(l) transformation 

($,$*) ^ (e^"$,e-^"$*), (2) 

where a is a constant phase. Since BEC takes place at very low energies, the interatomic 

potential can be approximated by the local potential f/(r — r') = 2X6{r — r') with A — ° 
related to the scattering length a. The Lagrangian simplifies to the form 

. l(r(M)%^-^*(M))-|lvr(M).V*(M) 

- V,,,(r)*-(r,t)<l>(r,t)-A(<l>*(r,t)<l.(r,())'. (3) 

The grand canonical system at the temperature T in the imaginary time formalism is defined 
by the Wick's rotated Lagrangian, 

L,m = r (r, r) (^-^ + - V^r) + ^(r, r) - A ($* (r, T)$(r, T)f , (4) 

where the chemical potential /i is introduced in order to consider states with an indefinite 
number of particles. Its value is to be adjusted so that the expectation value of the number 
operator is equal to JV, corresponding to a fixed particle density, n — j^- The Euclidean 
generating functional is defined as a path integral, 

Z[J] = J D^D^* g-i^f dr/d^r[L^[*]+J*(r,r)*(r,r)+J(r,r)**(r,r)] 

over the functions $(r,r) with a period P = in r. The partition function at thermal 
equilibrium is given by 2'[0]. The generating functional for connected Green's functions 
W[J] — \n.Z[J]. The effective action functional is given by the Legendre transform 

r[$] = W[J] - J dxJ{x)<^*{x) - J dxr{x)<l?{x), (6) 

where x — (r, r). The background field = — {^{x))j is a vacuum expectation of 
the quantum field operator in the presence of external source, J{x) — +ij2{x)). 
Performing the calculation we will use the two real components of the complex field = 
^ (01 (x) + i(f)2{x)) as independent variables. As a Legendre transform the effective action 
fulfils 

TTT^^-jii^) where i = l,2. (7) 

d(()i{x) 



4 



The physical value of the background field corresponding to J(x) = is thus 

determined by the stationarity equations 

— — - =0 where i = l,2. (8) 

Green's functions are obtained by functional differentiation of generating functionals. The 
one-particle Green's function (propagator) matrix reads 



Gij{x,y) 



6^Z 



6^W 



j=o 



Z[0] 6Ux)6j,iy) 

By differentiation of the effective action the one-particle irreducible Green's functions (proper 
vertices) are generated. The proper vertex 

-1 



6^T 








6(f)i{x)(j)j{y) 


$=$(0) 




J=0- 



G-Mx,y). 



(9) 



One-particle excitations, related to the poles of the full propagator, are determined by zero 
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modes of the Fourier transform Tij{p) = J dxe ^\(f,-(x)(j) {y) 
{uj,p). The matrix can be written in the form: 



<I>=<I>(0) 



where p stands for 



r(p) 



/ pi 

2m 



V 



11*2) with 



- yU + IIll UJ + Uu 

^ + n2i |^-^ + n2 

where Hn = Sn + i (S12+ 2*2), II22 = Sn - i (S12+ ^2) and Uu = | (S 
Ell = S22 being the normal self-energy and S12 = the anomalous one. 

In the following we study a homogeneous system, Vea;t(r) 0, when the background field 
$*^°''(x) = In this case the effective potential 

^\-^]\^{x)=^=const 



^^^^ PJd^r 
is an useful tool, since can be determined by the stationarity equation 



dV 



dV 



${0) 



0. 



(10) 



(11) 



$(0) 



The H-P theorem can be easily demonstrated in the effective action approach. Because of 
the invariance of L^[(l)i, ^2] under the transformation (j21), whose infinitesimal version is given 
by 



(12) 
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the generating functional [71,^2] is invariant under 





This imphes 



dx 



'SW 



SW 



-32 



-Ji 



= 0. 



(13) 



(14) 



Taking derivative over 02 (y) at "73^^!°^' '^^ which fulfils the symmetry breaking condi- 



tion 



<5r 



S4>i (x) 



$(0) 



0, we obtain 



/ 



dx y—(f)i{x) 
) / dx 



6^T 



502 (a;) 02(1/) 



+ 



6T 



3.(0) 



+ 02 (a;) - 



$(0) 



6(j)2{x)(j)2{y) 



S(pi{x) 

-0f )r22(^ = 0) = 0. 



{x)6(j)2{y) 



$(0)/ 



(15) 



$(0) 



Since 0f -* ^ 0, this means that r^^(a; = 0, = 0) = 0, there is therefore a zero- frequency 
excitation. The H-P theorem can be expressed in the form 

d^V 



= 0. 



(16) 



$(0) 



It is easy to observe that the gaplcss spectrum is just a consequence of the fact that the 
effective potential is a function of |$p and doesn't depend on 0i and 02 separately. This 
provides a useful criterium for the consistency of an approximation with the H-P theorem. 

For interacting atoms (a 7^ 0) the effective cannot be calculated exactly, so one resorts 
to approximations. The advantage of formulating the approximation for the effective action 
hes in the fact that all the Green's functions are obtained in a consistent way, through 
functional differentiation. 



III. THE LOOP EXPANSION 



The loop expansion is generated by calculating the path integral for Z[J] by the steepest 
descent method. Upon Legendre transformation the effective action is obtained as a series 
in h, whose power indicates the number of loops. The 1-loop effective action 

902 1 V72 1 V72 



A 



^0l(^+/.)01-^02(^+/.)02 



+ 



(01 + 02)" 



+ ^TrLnM[^], 



(17) 
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where M[$] : 



- + 3A0?(x) + \<Pl{x) + 2XMx)M^) 

-i^ + 2A0i(x)02(a;) -£ - + + 3A02(x) 



gives the Behaev-Popov approximation at finite temperature. Setting $(x) = $ yields 
the Toyoda's ^ result for 1-loop finite temperature effective potential 



V 



l~loop 



-/i|$|^+A|<l>|^ + Mi(^, |<l>| 



where 



(2vr) 



(18) 



(19) 



2- ■ P 

and ujk = ^ ~ + 4A|$p^ — 4A^|$|^. The H-P theorem is respected in the 1-loop ap- 
proximation, since the effective potential depends only on The false statement in 
Ref. 01 was due to wrong interpretation of excitation energy. One can verify explicitly that 



${0) 

^yl-loop 



0, if ^^^^ = -^{(f)f'\ 0) is a solution to the stationarity equation 



$(0) 



.^ + 2A|<l>(°f + ;iAJo(/i,|$(°f) 



l{0)_ 



-0. 



where 



d^k ^-/i + 3A|$|2 



2m 



(27r)3 

and the Bose-Einstein distribution function nsiuj) 



'l + 2nBiuJk)) 



(20) 
(21) 



(e,fl"_l) • 



The chemical potential /i can be eliminated in favor of the particle number density, using 
the relation 



n 



dV 
dfi 



|$oP + ^/o(/i, |$o| 



$(0) 



where 



^o(yU, |$| 



d^k 



2m 



-/i + 4A|<l>| 



(2vr) 



l + 2ni3(c^fc)) 



(22) 



(23) 



Renormalization can be carried out by dimensional regularization. Denoting the particle 
density in condensed state |$oP by uq and observing that in the leading order rig = ra, the 
particle density to the order h can be written as 



n = no + 



d^k £ + 2An 



(2vr)^ 



Bog 



l + 2n^(c.f^^ 



(24) 



where Bogoliubov frequency uj^°^ = \/^(^ + 4An). This expression has been used |3| to 
calculate the condensate depletion 

8(an)3/2 



uq = n 



3A 



m 
.2^, 



2/3 



C(3/2) - 27ri 



'2(3an 



m 



(25) 
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The depletion and other physical quantities (energy density, pressure, . . .) derived from 1- 
loop effective potential are in agreement with the Lee and Yang results to the lowest order 
in arv^/'^. However, the 1-loop results are dubious in the vicinity of the phase transition, since 
the higher-loop contributions are significant at such a temperature. 

IV. THE OPTIMIZED EXPANSION 

The optimized expansion consists in introducing an arbitrary parameter VL into the 
Lagrangian density 

L^i^, n] = r(r, r) [-^ + |-V^ + n\ $(r, r) 



+e 



<l>*(r,r)(/i-l])$(r,r)-A(<l>*(r,r)$(r,r))' . (26) 



For e = 1 the dependence on Vt cancels and the modified Lagrangian coincides with the 
original one (@]). Calculating Z[J] and performing the Legendre transform yields the effective 
action as a series in a formal parameter e. The exact result would not depend on however 
such a dependence appears in the n-th order truncation, r*^"^[$,f2], obtained after setting 
e = 1. We exploit this freedom by choosing the value of fi, which fulfils the minimal 
sensitivity requirement 

0, (27) 



5^ n°p' 

in the given order approximation r*^")[$, 1]°^*]. The optimal value of VL changes from order 
to order, improving the convergence properties of the scheme. The approach is equivalent 
to a systematic re-summation of the perturbation series. 

Gaussian approximation to the effective action is obtained in the first order of the opti- 
mized expansion. This yields the Gaussian effective potential 

y«[<l>, fi] = +A|<I>|^ + /i(fi, |<l>n + (/i - fi)/o(l^, |$|') + 3A/2(fi, |<I>H (28) 

with self-consistency condition 

f?i^ = |(,-n)+6A/„(a|4ra§ = 0. (29) 

where Ii and Jq are defined by (fTT?|) and ipHjl . respectively. Since both the effective potential 
and the optimization condition depend only on one clearly sees that the H-P theorem 
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is respected. The symmetry is spontaneously broken at |$p = no which renders the effective 
potential stationary: 



= -fi+2Xnl + XJo{Q,no) + [{fi - Q) + 6XIo{Q,no)] 

1*1^0 











2 



d\^ 

with Jo given by (j25)- The particle density in our approximation reads 

dVW 



= (30) 

1*1 



n 



no + Io{^,no) (31) 



dfi 

Upon eliminating /i by (j^ . the EqslHUl and IMl simplifv to 

n = no + Io{n,no) 

Q = 2Xnl + AJo(fi, ^o) + 6A/o(l^, no) (32) 
which, after numerically eliminating Q, determine no{n,i3). 

V. RESULTS AND CONCLUSIONS 

The critical temperature for BEC, Tc, is the temperature below which the symmetry is 
spontaneously broken. Its value can be calculated from the condition that |$p = no = at 

2 

the phase transition. For the ideal gas T^'^ = (^) ( ) , but calculating the shift of the 



critical temperature in the presence of interactions 



= (33) 

c 

generates controversy, even for the dilute gas. Different powers of the leading behavior 



in an^^^ with different coefficients have been reported by various authors ij]. is usu- 
ally derived from the effective 3-dimensional theory, arguing that only the zero Matsub- 
ara modes determine the critical behavior. In this approximation the leading behavior 
ATc ~ 1.32an^/^ and next-to- leading order corrections have been determined 

Here we show the approximate results for Tc derived in the original theory in (3 + 1)- 
dimensions. The 1-loop approximation for Tc has been calculated [16| by setting the value 
of the background field Hq to zero in Eq. |^ which yields the equation 

n = Io{2Xn,n), (34) 
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to be solved numerically. We calculate the Gaussian approximation for by setting no = 
in Eq. EH which results in the coupled pair of equations 



n = AJo(fi,0) + 6A/o(fi,0). 



(35) 



In Fig.l the numerical results for are compared with the results derived in the dimen- 



sionally reduced theory 
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, 3] • For the experimental evaluation of the accuracy achieved 
in homogeneous systems is still insufficient. The theoretical results differ much, even at small 
values of an^^^. It is interesting to observe that the leading behavior in 1-loop approximation 
is of the form AT, 

4v^ VanV3 



3 [C(3/2)]2/3 

while the Gaussian approximation gives AT, 



Vj; Van^ /j^ ^^j^]-^ j^Y^Q gjgj^ opposite to that obtained by Toyoda 



with the coefficient -\/2 times 



3 [C{3/2)]2/3 

larger than 1-loop result. In both cases the square root dependence on an^^^ is obtained, 
similar to Lee and Yang result 0] , but in difference with the linear dependence derived from 
the dimensionally reduced theory. The behavior of in the Gaussian approximation at 



Tom the 1-loop behavior and there is 
"1 

16| by the 1-loop approximation. It 



larger values of an^^^ becomes qualitatively different 
no sign of the reentrant phase transition suggested 
would be interesting to investigate this issue in higher order approximations. Both the loop 
expansion, as well as the optimized one, offer a systematic and consistent way for such a 
study, without violating the H-P theorem. 




0.8 ,1 



FIG. 1: The critical t emp erature for BEG in the Gaussian approximation(so/id line), compared 
with the 1-loop result {dotted line), and the results of the 3-dimensional theory: the leading 
order in an^/^ j3l {dashed line), and next-to leading order Q| {dotted- dashed line). 
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